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We study the sensitivity of the microwave and submillimeter transitions of the isotopologues of 
hydronium to the variation of the electron-to-proton mass ratio (i. These sensitivities are enhanced 
for the low frequency mixed inversion-rotational transitions. The lowest frequency transition (6.6 
GHz) takes place for isotopologue H2DO + and respective sensitivity to ^-variation is close to 200. 
This is about two orders of magnitude larger than the sensitivity of the inversion transition in 
ammonia, which is currently used for the search of /^-variation in astrophysics. 
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I. INTRODUCTION. 

The search for the possible variation of fundamental 
constants is a poor man's way to look for physics beyond 
the Standard Model. There are two main approaches 
to such studies. First is to make laboratory measure- 
ments with different atomic and molecular clocks in or- 
der to look for local variations on the time scale of the 
order of a year. Second is to compare laboratory mea- 
surements with astrophysical observations, where one can 
study variations on a much bigger space and time scale 
up to 10 billion years. Advantage of the first approach 
is in unprecedented accuracy of experiments with atomic 
clocks. On the other hand, the second approach allows 
exploring bigger scales and different environments. For 
example, the local matter density in the laboratory ex- 
periments and in interstellar molecular clouds differ by 
more than 10 orders of magnitude. We can conclude that 
these two approaches are complementary to each other 
In this paper we do not discuss laboratory experi- 
ments and focus on astrophysical applications. 

In astrophysics all frequency shifts with respect to 
the laboratory values are conventionally interpreted as 
Doppler shifts, Au/uj — V/c, where V is the line of sight 
velocity of the object and c is the speed of light. How- 
ever, if we observe several lines from the same object, 
we can study (small) velocity offsets between them. In 
the search for the variation of fundamental constants we 
look for the correlation between these offsets and sen- 
sitivities of the observed transitions to the variation to 
fundamental constants. 

It is now well known that the tunneling transition 
in ammonia is highly sensitive to the electron-to-proton 
mass ratio fj, = m e /m p [2J. This transition is often 
observed from the interstellar dark molecular clouds in 
our Galaxy. Comparison of the apparent Doppler shift 
for this transition with those of rotational transitions of 
other molecules provides a very sensitive test for the de- 
pendence of the mass ratio /i on the local matter density 



Moreover, ammonia was detected in two high red- 
shift objects with z = 0.68 and z = 0.89. That allows to 
place the most stringent present limits on /i- variation on 
a cosmological time scale [(HU . 

Because all observed ammonia transitions have the 
same sensitivity to the variation of the fundamental con- 
stants, one has to use rotational transitions in other 
molecules as a reference. This can lead to the system- 
atic effects caused by different spacial distribution of the 
molecules in gas clouds (so called Doppler noise) [Toj . 
Therefore, it is desirable to find molecules, where one can 
simultaneously observe several lines with high and differ- 
ent sensitivities to /i-variation. In this case the Doppler 
noise is minimized because all lines are observed from the 
same gas. By now we know several such molecules, in- 
cluding OH [ll|, 0j| and CH [l3| , partly deuterated am- 
monia [l4l| . h y dronium ion (H30 + ) and methanol 
(CH 3 OH) 111 El. These molecules often have higher 
sensitivity than ammonia and can be used as indepen- 
dent source of information on variations of fundamental 
constants. 



In this paper we consider the hydronium isotopologues 
H 3 0+, H 2 DO+, HD 2 0+, and D 3 0+. Like ammonia, 
hydronium also has a double minimum vibrational po- 
tential. The inversion transitions occur when the oxy- 
gen nucleus tunnels through the plane of the hydro- 
gen/deuterium nuclei. This leads to an inversion split- 
ting of the rotational levels. The splitting in HaO + is 
very large, about 55 cm -1 , as compared to the 0.8 cm -1 
splitting in NH3. Because of that the inversion transition 
is effectively mixed with rotational transitions and the 
inversion-rotational spectrum of hydronium is observed 
in the submillimeter-wave region. Here we show that 
all hydronium isotopologues have mixed transitions with 
high sensitivity to /i-variation and we also give improved 
estimate of the sensitivity coefficients for the main iso- 
topologue H 3 + . 
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II. GENERAL FORMALISM 
A. Sensitivity coefficients 

The dimensionless sensitivity coefficients to the varia- 
tion of fundamental constants can be defined as 14] 



Suj 5a 5[i Sg n 

u> a [i g n 



(1) 



Here a « 1/137.036 is the fine-structure constant, /i = 
m e /m p ~ 1/1836.15, g n is the nuclear (/-factor, and Q a , 
Qn, and Q g are the corresponding sensitivity coefficients. 
Note that g n is not a fundamental constant, but it weakly 
depends on the quark masses [l8[ and has to be consid- 
ered as independent parameter whenever magnetic hy- 
perfine structure is involved. 

It is known that tunneling and rotational transitions in 
molecules, built from light elements, are mostly sensitive 
to fi- variation: > 1, while Q a <C 1 and Q g <C 1. 
In the following we disregard coefficients Q a and Q g 
and concentrate on calculating the dominant sensitiv- 
ity coefficient Q p . In the next two sections we consider 
purely inversion transitions. Note, that for the asym- 
metric isotopologues H2DO + and HD20 + such transi- 
tions are not observable; they occur only in combination 
with rotational transitions w r [l4| . We will discuss mixed 
inversion-rotational transitions in Sec. Ill Dl 



B. Inversion potential and reduced mass 

As explained above, we can link the variation of the in- 
version transition (5u;j nv /cL>i nv to the variation of the single 
fundamental constant /x: 



■ ^ 



(2) 



We can present the Hamiltonian for the inversion process 
as (if not said otherwise, we use atomic units throughout 
the paper): 



(3) 



where x is the tunneling coordinate and M is the respec- 
tive reduced mass. For example, x can be the distance 
from oxygen to the plane of hydrogens. Following (l9j we 
parameterize the potential U(x) as follows 



U(x) = -kx 2 + -dx* 
w 2 2 



(4) 



Parameters k, 6, c, and d can be found by fitting vibra- 
tional frequencies for the molecule H 3 + . As a result 
we obtained k = 0.08917, b = 0.04977, c = 1.36954, and 
d= 0.00138. 

Another way to find the parameters in Eq. (QJ is to fit 
the inversion frequency of the ground vibrational state 



and the barrier height. The latter was accurately deter- 
mined within "semi-empirical" approach in Ref. |20j to be 
652.9 cm- 1 . In this case we got k = 0.08907, b = 0.04967, 
c = 1.30429, and d = -0.00541. To find the barrier 
height AU = U max — U m - ln we need to solve equation 
dU(x)/dx = 0. Obviously, C/ max = U(0) = b = 10901 
cm . As to :c m i n and U m - ln , we need to solve numeri- 
cally the equation 

k + 2d a4 in -26 c e~ ca; ™ = , (5) 

which gives x min — 0.5635 and f7 m i n = 10248 cm -1 . 
From this we get AU « 653 cm -1 in a perfect agreement 
with the semi-empirical result 652.9 cm -1 from Ref. [20l ] . 

Comparing two sets of parameters we see that they 
are very close with the exception of the parameter d. We 
have checked that both sets of parameters for the poten- 
tial U (x) lead to very close results for the sensitivity co- 
efficients. The numerical results which will be discussed 
in the following are obtained with the second set. 

In order to find eigenvalues of the Hamiltonian (J3]) we 
need to know the reduced masses for different isotopo- 
logues of hydronium. For the symmetric species HsO + 
and D3O" 1 " three protons (or deuterons) can be treated 
as a single particle with the mass m = 3 m p ^ . Then the 
reduced mass for the inversion mode can be found within 
two simple models. 

If we assume that the H-H distances are constant dur- 
ing inversion, then the vibration involves only a change in 
the distance between the plane of hydrogens and oxygen. 
The reduced mass M(H30 + ) can be written as 



M(H 3 0+) = 



3m p Mo 
3m p + Mq 



(6) 



where Mo is the mass of the oxygen nucleus. 

Alternatively, we can assume that the H-0 distances 
are constant. Then the inversion coordinate corresponds 
to the angle and the reduced mass is given by [19j 



M(H 3 0+) = 



3m p (Mq + 3m p sin 2 9) 



3m„ + Mi 



o 



(7) 



Here 9 is the angle between the plane of the hydrogen 
atoms and an H-0 bond. Using the value of the angle, 
a, between two H-0 bonds a(HOH) = 111.3° 13 we 
obtain 9 w 17.6°. To find the reduced mass of the DsO + 
molecule we only need to change m p to rrid 2m p in 
Eqs. ([6]) and (JT]). Considered here models present two 
limiting cases. The difference between them constitutes 
only 1.7% for H 3 0+ and 3.4% for D 3 0+. More accurate 
values for reduced mass should lie between these limits. 

The motion of the asymmetric species H2DO" 1 " and 
HD2 + is more complex and finding reduced masses is 
more difficult. We can consider it as a free parameter 
that can be found from fitting the theoretical inversion 
frequency of the ground state to the experimental value. 
Note that the inversion frequencies of the ground states 
are measured with high precision for all the molecules 
H 3 0+, H 2 DO+, HD 2 0+, and D 3 0+ 0. 
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TABLE I: The reduced masses M, the experimental and the- 
oretical inversion frequencies ajj nv , and the numerically calcu- 
lated sensitivity coefficients Q M ,i n v 







(•)■,„,. (cTCl' 
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Molecule 


lvl ^a.u.J 


Experim. 


i iieor . 


d 


H 3 0+ 


4639 


55.35 


55.79 


1.97 




4657 




55.35 


1.98 




4718 




53.94 


1.99 


H 2 DO+ 


5419 


40.52 


40.52 


2.15 


HD 2 + 


6485 


27.03 


27.03 


2.36 


D 3 0+ 


8012 


15.36 


15.92 


2.65 




8120 




15.36 


2.67 




8287 




14.54 


2.70 



TABLE II: The reduced masses M, the ground state vibra- 
tional energy Eo, the action over classically forbidden region 
5, and the sensitivity coefficients Q Ml mv obtained in the WKB 
approximation. 



Molecule 


M (a.u.) 


Eo (cm" 1 ) 


S 


Q fi.inv 


H s O + 


4639 


339 


1.36 


1.92 




4657 


339 


1.36 


1.91 




4718 


337 


1.36 


1.90 


H 2 DO + 


5419 


322 


1.62 


2.10 


HD 2 0+ 


6485 


302 


1.96 


2.33 


D 3 0+ 


8012 


279 


2.45 


2.64 




8120 


278 


2.44 


2.63 




8287 


276 


2.44 


2.61 



C. Sensitivity coefficient for inversion transition 

A numerical integration of the Schrodinger equation 
with the potential ((4]) allows us to find eigenvalues of the 
Hamiltonian and, respectively, to calculate the inversion 
frequency. Then we can find the sensitivity coefficients 
Qfj,,mv by numerical differentiation taking into account 
that reduced mass M in atomic units scales as (J, . 

In Table UJ we present the results of the numerical cal- 
culations for the ground state inversion frequencies and 
the sensitivity coefficients for the isotopologues of hydro- 
nium. For the symmetric species H3O" 1 " and DsO + we 
present results for three values of the reduced mass. The 
first and the third values correspond to Eqs. (HJ) and (J7J, 
correspondingly. The second value is found by fitting the 
ground state tunneling frequency. As seen from the ta- 
ble the coefficients Q^^m, are almost insensitive to these 
small changes of the reduced masses. Even for the heav- 
ier isotopologue D 3 + the differences between calculated 
values of Q^inv do not exceed 2%. 

Another way to find the sensitivity coefficients Q^inv 
is based on the semi-classical Wentzel-Kramers-Brillouin 
(WKB) approximation. A detailed description of this 
approach and its application to calculations of Q^inv for 
different molecules was repeatedly discussed earlier (see, 
e -g-, [Still Gil [III). For this reason we present here only 
the main expressions. 

Following Landau and Lifshitz [23| we can write the 
inversion frequency as 



TABLE III: The recommended values of the sensitivity coef- 
ficients Q Ml inv The uncertainties are given in parentheses. 



2E, 



°e- s , 



(8) 



where S is the action over classically forbidden region 
and E is the ground state vibrational energy calculated 
from the bottom of the well U m i n . Harmonic approxima- 
tion 2Eq = oj v (uj v is experimentally observed vibrational 
frequency) is not applicable to hydronium and Eq has to 
be found by solving Eq. (j3|) numerically. After that using 
the experimentally known a>i nv we can obtain the action 
S from Eq. ©. 

It was shown in [(| [l5j] that the sensitivity coefficient 





H 3 0+ 


H 2 DO+ 


HD 2 0+ 


D 3 0+ 




2.0(1) 


2.2(1) 


2.4(1) 


2.7(1) 



?M,inv can be expressed through S, E , and AU as follows 
S SE 



Q 



1 



fi , mv 



2(AU - E ) 



(9) 



The third term in Eq. (j9|) was first obtained in Ref. [f| ■ 
It was shown there that the numerical factor for this term 
depends on the form of the potential barrier. (It is worth 
mentioning that there is a typo in Eq.(12) of @: the 
coefficient 1/4 should read as 1/2.) The square and tri- 
angular potential barrier were considered. For both of 
them the coefficients in front of the third term were cal- 
culated and the average was taken. This average value is 
reproduced in Eq. (|9]). Here we considered a more real- 
istic, parabolic form of the potential barrier and found, 
that solution for this case exactly coincides with Eq. (JSJ . 

Knowing the values of S, Eo, and AU we can calcu- 
late the sensitivity coefficients Q Mi i nv in the WKB ap- 
proximation. The results are presented for all molecules 
in Table HU The WKB approximation is known to work 
well for S > 1. As seen from Table HU for hydronium 
S is close to unity. For this reason we cannot anticipate 
high accuracy for the sensitivity coefficients obtained in 
this approximation. Nevertheless, comparing the results 
of Tables Q] and [Til we see a reasonable agreement be- 
tween the numerical and semi-classical values of Q^inv 
The largest discrepancy does not exceed 5%. Based on 
the comparison of all approximations, considered here, 
we present the recommended values of the sensitivity co- 
efficients Q Mi inv m Table Mil We estimate the accuracy 
of our calculations to be not worse than 5%. 

The coefficient Q M ,i n v for H 3 + calculated in this work 
differs by approximately 20% from that in [l^]. In 
Ref. [HI this coefficient was found in two ways: in the 
semi-classical approximation and from the plot where the 
inversion frequency was presented as a function of the 
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reduced mass of hydronium isotopologues. The reduced 
masses were estimated on the basis of Ref . [2(| ■ 

We already mentioned above that the accuracy of the 
semi-classical approach for H3 + is not very high as the 
action S ~ 1. In addition, the answer strongly depends 
on the value of Eq in Eq. (|9|). Here wc found Eg solv- 
ing an eigenvalue problem for Eq. @, while in Ref. [15] 
Eq was extracted from Ref. [20J to be approximately 
equal to 400 cm -1 . The second, graphical method used 
in [r| strongly depends on the values of the reduced 
masses, which were approximately determined from Ref. 
[20| . These approximations caused overestimation of the 
Q-factor. 

Here, in order to improve the accuracy, we have done 
numerical calculations within the simple one-dimensional 
model described above. Optimizing parameters of the po- 
tential U (x) we were able to reproduce the experimental 
low-lying energy levels (the first set of the parameters) 
and the transition frequency and the barrier height (the 
second set of the parameters) with an accuracy better 
than 1%. However, we still estimate the accuracy of this 
simple model for the Q-factors to be about 5%. 



D. Sensitivity coefficients for mixed transitions 

Sensitivity coefficients for mixed inversion-rotational 
transitions in H30 + were considered in [Fa ]. Here we 
extend this discussion to all four isotopologues of hydro- 
nium. Rotational motion is described by the Hamiltonian 
of the non-rigid asymmetric top. We write it in the ba- 
sis set \J,K,s) of the oblate symmetric top (K = Kc)- 
Additional quantum number s corresponds to the sym- 
metric (s = +) and asymmetric (s = — ) inversion state. 
Rotational parameters A s , B s , and C s depend on the 
quantum number s. This dependence contributes to the 
inversion transition + <H> — and can be considered as 
centrifugal corrections to the inversion frequency F. 

The effective inversion-rotational Hamiltonian is diag- 
onal in quantum numbers J and s and mixes states with 
K and K ± 2. Matrix elements diagonal in K have the 
form: 



H 



K.K 



-F 



(10a) 



+ - {A s +B S ) [J(J+ 1) -K 2 ] + C S K 2 (10b) 

+ Aj [J (J + l)] 2 - A KJ J(J + l)K 2 - A K K A , (10c) 

where the last three terms describe centrifugal correc- 
tions to rigid rotor. Matrix elements nondiagonal in K 
are given by: 



H k ,k+2 = \{A s - B s ) [J(J + 1) - K(K + 1)] 
x [J(J+1)-(K +1)(K + 2)} 1/2 . 



1/2 



(lOd) 



Effective Hamiltonian (|10[) has ten parameters, which 
have to be fitted to the experimental spectrum. For 



the symmetric species we require that A s — B s and 
A ^, / = 0, leaving seven free parameters. Such fits were 
done for H 3 0+ in Ref. [H, for H 2 DO+ in 0,111 HI, 
for HD 2 0+ in ^Hl], and for D 3 0+ in In the cited 
literature the centrifugal terms (I10c|) are also assumed to 
depend on s. These s-dependent corrections appear to 
be much smaller than all other parameters of the effec- 
tive Hamiltonian and here we reduce the number of free 
parameters by neglecting them. 

In order to calculate sensitivity coefficients for mixed 
transitions we need to specify how parameters of the ef- 
fective Hamiltonian depend on ji. The inversion transi- 
tion energy is given by the parameter F, whose depen- 
dence on /j, was discussed in Sec. Ill CI Averaged rota- 
tional constants A, B, and C depend on the respective 
moments of inertia and are, therefore, proportional to /x 
(A = (A + + A~)/2, etc.). Centrifugal parameters Aj ap- 
pear in the second order of adiabatic perturbation theory 
and are quadratic in fi. 

It is more difficult to determine the /^-dependence of 
the differences A + - A~ , B+ - B~ , and C + - C~ , which 
present centrifugal corrections to the inversion transition 
energy F. All of them must have the same dependence 
on n, so we will consider only Wjc = C + — C~ . 

Clearly, the rotational constant C generally depends 
on the inversion coordinate: C = Co + C\x 2 + . . . . 
This generates a correction to the inversion potential ((H) , 
U(x) — > U(x) + C\x 2 K 2 . To a first approximation we 
can substitute x with its equilibrium value a; m i n . Con- 
sequently, C/ m in — > U m i n + Cix^^K 2 . We can now use 
Eq. (jSJ| to estimate the respective change in inversion 
frequency Wi nv — > w; nv + WkK 2 . After some algebra we 
arrive at: 

S 

Wk = Winv AU-E ClX ™™- 
Differentiating this expression in respect to fi we get: 



sw. 



K 



w. 



K 



= Q 



1\ Sji 



(12) 



We conclude that centrifugal corrections to inversion fre- 
quency scale as jjQ^,^+ 1 / 2 1 which is sufficiently close to 
the estimate [6|. 

Since we determined how the parameters of the effec- 
tive Hamiltonian depended on /i, we can find the sen- 
sitivity of the mixed transitions by numerical differen- 
tiation of the eigenvalues of the effective Hamiltonian. 
Results of these calculations are listed in Tables IIVI and 
fVl Note that if we neglect all centrifugal corrections, we 
can write the approximate expressions for the frequency 
and the sensitivity of the mixed inversion-rotational tran- 
sition [l3 |: 

(13a) 
(13b) 



k^mix LUr i ^inv , 



k^mix ^mix 

For small values of the quantum numbers J and K con- 
sidered here the centrifugal corrections are small and ap- 
proximation (|13p is rather good. 
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TABLE IV: Sensitivities of the low frequency mixed inversion- 
rotational transitions in symmetric isotopologues of hydro- 
nium. Molecular states are labeled with quantum numbers 
Jk. Error bars for sensitivity coefficients Q M correspond to 
the errors for inversion transitions in Table [Till 



Transition 


Frequency (MHz) 




Upper 


Lower 


Theory 


Exper. 








T T C\ J r 






i — 
L i 


9+ 


QA7A79 

OU f U / Z 


QA71 QO /I 
OU / iyz.4 


6.4(5) 


■'>+ 
6 2 


2 2 


OOOU40 


O04/ y / .4 


-3.5(5) 


■'>.+ 


2i 


QQQ1 en 


ooo4oo.o 


-3.1(4) 


q+ 
d o 


9 — 
2() 


oy / iyo 


oyoz / z.4 


-3.0(4) 


u o 


1 + 
->-() 


yo4oyU 


yo4/ 11. y 


2.7(2) 


A+ 


O — 

3 3 


1 aoi ftft/l 
lUol004 


lUol2yo. ( 


-0.6(2) 


4 2 


q 

J 2 


1 A71 1 C/i 
1U / 1 104 


iuoyozo.0 


-0.5(2) 


J 2 


d 2 


1 (?oi QOfi 
lOZloZO 


idzi / oy.u 


2.0(1) 




9+ 






2.0(1) 


i-r 


If 


1655832 
D 3 + 


1655833.9 


2.0(1) 




377 


450608 


450709.7 


2.7(1) 






454910 


454940.3 


2.7(1) 


2-r 


n 


456194 


456211.4 


2.7(1) 


i-r 


it 


459918 


459917.7 


2.7(1) 


22" 


■21 


460496 


460493.1 


2.7(1) 


3 3 " 




462080 


462074.6 


2.7(1) 






120117 




7.5(4) 


2i 




220408 




-2.5(2) 


2,! 


1,7 


221782 




-2.5(2) 



III. DISCUSSION 

Our final results for sensitivity coefficients for 
inversion-rotational transitions between low-lying states 
of all hydronium isotopologues are listed in Tables llVl and 
IVl There we also give calculated frequencies and experi- 
mental frequencies, where known (see Refs. (25l. [28l - [30| ). 
The errors for the Q-factors correspond to the uncertain- 
ties in Table Hill For low- frequency mixed transitions the 
errors are increased due to cancelations between rota- 
tional and tunneling contributions in Eq. (|13b|) . For the 
experimentally observed transitions we use experimental 
frequencies to calculate Q-factors. For unknown transi- 
tions we have to use predicted values. The accuracy of 
our model with limited number of centrifugal corrections 
decrease with increasing rotational energy. For given J 
the latter is smaller for heaviest isotopologue. Therefore, 
the accuracy of the model is highest for D 3 + (~ 100 
MHz) and lowest for H.30 + , where the maximum error 
is almost 1 GHz. However, the errors in Q-factors from 
the inaccuracy in frequency do not exceed 2% and can be 
neglected compared to the error associated with Q^inv 

It is clear from Eq. (|13j) that mixed transitions can 
have high sensitivities only if |w m i x | <C Wj nv [note, that 
(|13al) allows negative frequencies w m ix]- Thus, we are pri- 
marily interested in transitions from the microwave range 
and low-frequency submillimeter transitions. One can see 
from Tables IIVI and [V] that all isotopologues have such 



TABLE V: Sensitivities of the low frequency mixed inversion- 
rotational transitions in asymmetric isotopologues of hydro- 
nium. Molecular states are labeled with quantum numbers 
T 3 



Transition 


Frequency (MHz) 




Upper 


Lower 


Theory 


Exper. 








H 2 DO+ 






3+ 

°1,2 


2(3,2 


210994 


211108.8 


-5.9(6 


io.i 


277,1 


250920 


250914.1 


6.8(5 


0+ 

J 2,l 


2 1,1 


611161 


010001 
012831.8 


—3.6(4 


0+ 
^2,2 


2l,2 


412156 


A 1 <~» 1 1A O 

412130.2 


—2.5(3 


0+ 
°3,0 


z 2,0 


632 


oo2yui. ( 


—1.2(2 


n + 


22,1 


649742 


649653.4 


— 1.2(2 


a - 

°0,0 


1 + 

1 i,o 


(i7QOOn 


( 6201 .U 


3.2(2 


A + 
4 2,2 


°1,2 


/ lo9oo 


/ 1082/.9 


1 A/1 

— 1.0(2 


A + 

4 i,3 


3o,3 


716961 


171 PAr A A 

716959.4 


— 1.0(2 




9+ 
z 2,l 


6633 




219(18 


\o 


9+ 


51108 
HD 2 0+ 




29(2 


0o7o 


K 


380753 


380538.0 


4.0(2 


0+ 

°2,2 


2l,2 


474792 


A T A f A 1 1 

474541.1 


-1.4(2 


^2,3 


5+ 

°3,3 


525419 


525451.5 


3.1(2 


3+ 
°3,0 


22,0 


610712 


610573.1 


-0.8(1 


3+ 

°3,1 


22,1 


626907 


627069.8 


-0.8(2 


6 2,4 


6j 4 


632532 


632776.9 


1.9(1 


3l,2 


3+ 

°2,2 


634305 


633793.2 


2.8(1 


4 r, 3 


4,3 


707925 


707552.6 


2.6(1 


!o7i 




728086 


728420.2 


2.6(1 




Ki 


35047 




33(2 


2+ 

z 2,0 




93734 




-11(1 


2+ 

z 2,l 


17.1 


129633 




-7.7(6 



transitions with both positive and negative sensitivities, 
whose absolute values are significantly larger than unity. 

The lowest frequency 6.6 GHz corresponds to the tran- 
sition 17" j 2%^ in H 2 DO+. This frequency is 184 
times smaller than respective inversion frequency. Con- 
sequently, the sensitivity coefficient for this transition is 
close to two hundred. Other frequencies are at least an 
order of magnitude larger and respective sensitivities are 
much smaller. 

As we pointed out above, our present value of Q^inv 
for H 3 0+ is 20% smaller than that in Ref. [Hj]. Conse- 
quently, the absolute values of our Q-factors for mixed 
transitions are also significantly smaller. Still, the dif- 
ference AQ between sensitivities of the 307 GHz line on 
the one hand and 365 GHz and 396 GHz lines on the 
other hand is close to 10, or 3 times larger, than AQ 
for ammonia method @. These three lines have been 
observed from the interstellar medium (see (3l| and ref- 
erences therein) . That makes hydronium very promising 
candidate for the //-variation search. 

Submillimeter spectra of H 3 + are usually observed 
from warm and dense star forming regions, where hydro- 
nium is one of the most abundant molecule [32(. Up to 
now deuterated hydronium was not observed from the 
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interstellar medium. However, deuterated isotopologues 
of ammonia were observed several times [33|. Thus, it 
is possible that hydronium isotopologues can be found 
in the future. As we have showed here, this can give 
additional possibilities to study /i-variation. 
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